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Abstract 

We study domain structures in bistable systems such as the Ginzburg- 
Landau equation. The size of domains can be controlled by a global 
negative feedback. The domain-size control is applied for a localized 
spiral pattern. 
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Spatially localized states have been observed experimentally in binary 
fluid mixtures [Q], electroconvection in nematic liquid crystals j|] and in 
granular media undergoing the Faraday instability ||. Some simple model 
equations have been studied to understand the mechanism of the localized 
states found in dissipative systems. Soliton-like localized states have been 
found in the quintic complex Ginzburg-Landau equation and the coupled 
complex Ginzburg-Landau equations [|]. Wormlike localized states were 
studied with the anisotropic complex Ginzburg-Landau equation coupled 
with a scalar mode || . A self-trapping mechanism works for localized states 
in the quintic Swift-Hohenberg equation Q. The long-range inhibition is 
important for localized states in some reaction diffusion equations ||. On 
the other hand, controlling chaotic dynamics has been investigated with the 



OGY method and the feedback method ||, 10 1. Zykov et al. studied the 



control of spiral waves in a spatially extended system by a global feedback 
0- 

We study the control of the domain size of localized domains in spa- 
tially extended bistable systems. Our first model equation is based on the 
Ginzburg-Landau equation coupled with an inhibitory medium. 

du 3 d 2 u 

dt dx 2 ' 

dv d 2 v 
T _ = _„ + ir( „_ c) + D _, (1) 

where u obeys the Ginzburg-Landau equation, v denotes the inhibitory vari- 
able, D is the diffusion constant of v, r is the time constant of v and c is 
a parameter between -1 and 1. This model equation is a reaction diffusion 
equation. The system size is L and the periodic boundary condition or the 
no-flux boundary condition is assumed. If D is sufficiently large, a localized 
state is obtained as in [8]. If D is infinitely large, v is uniform for < x < L, 
that is, v = (v) = (l/L) J Q L vdx. The second equation in Eq. (1) is reduced 
to 

r^ = -(v)+K((u)-c), (2) 

where (u) = (l/L) J Q L udx. If the adiabatic approximation is assumed, (v) = 
K((u) — c). Substitution of this relation into the first equation of Eq. (1) 
yields 

du , ,,,, . . d 2 u , . 

_ =u _ u 3_ K((u) _ c) + _ (3) 

where the third term on the right-hand side represents the global negative 
feedback. If the third term —K((u) — c) takes a small constant value b, there 
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are two stable states near u ~ ±1. If the initial condition takes a domain 
structure, the domain wall moves with a constant velocity for nonzero b. 
The domains of positive (or negative) u become dominant if b is positive (or 
negative). In our negative feedback model with K > 0, the effective control 
parameter —K((u) — c) decreases (increases) as the domain size of positive 
(negative) u increases. Finally, (u) = c is attained and the domain growth 
stops. Then, the size of the domain of u = ±1 is approximately (1 ± c)L/2. 
Namely the parameter c determines the domain size. We can control the 
domain size by changing the parameter c. Figures 1(a) and (b) display the 
time evolution of the domain structure at K = 0.5, c = —0.2 and L = 200. 
The initial condition is u(x) = 1 for 98 < x < 102 and u = — 1 for x < 98 and 
x > 102. The numerical simulation was performed with the pseudospectral 
method of mode- number 1024 and timestep 0.01. Figure 1(a) displays the 
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Figure 1: (a) Time evolution of u(x,t) by Eq. (1) at r = 1, K = 0.5, c = 
-0.2, D = 10000 and L = 200. (b) Time evolution of u(x, t) by Eq. (3) at 
K = 0.5, c = -0.2 and L = 200. 



time evolution for Eq. (1) with D = 10000 and r = 1. A localized state with 
size 75.7 is obtained. The domain size is calculated as a width of the region 
of u > 0. Figure 1(b) displays the time evolution by the Ginzburg-Landau 
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equation (3) with the global negative feedback. A domain with a fixed size 
is obtained as a stationary state. The final size of the domain u = 1 is 80. 
Since the diffusion constant D = 10000 in Eq. (1) is very large, the time 
evolution by Eq. (1) is close to that by Eq. (3). The domain size is well 
approximated at the value (1 + c)L/2 = 80. 

We can control the domain size even if the system has three stable states. 
The model equation is the quintic Ginzburg-Landau equation: 



du 



-— = —au + u 



u° + e + 



u 



dx 2 ' 



(4) 



where a and e are parameters. If e = and a = 3/16, the potential energy 
U = au 2 /2 - n 4 /4 + n 6 /6 - eu takes the same local minimum value at 
u = and ±\/3/2. Two kinds of domain walls which connect — y/3/2 and 0, 
and and \/3/2 do not move at the parameters e = and a = 3/16. For the 
other parameter values, the domain walls have finite velocities. We assume 
a model with a global negative feedback as 



— = u 3 - u 5 - ^{(u) 



d 2 



u 



ci 



K 2 ((u 2 )-c 2 )u + -^, (5) 



where (u) = (1 / L) Jq udx , (u 2 ) = (1/L) Jq u 2 dx, and Ki,K 2 are coupling 
constants and c\ , c 2 are the control parameters which determine the domain 
sizes. Since the parameter value (a, e) = (0, 0) in Eq. (4) is a codimension 
2 point, we need two types of negative feedback. If the negative feedback 
succeeds, the final state satisfies (u) = c\ and K 2 {{u 2 ) — C2) = 3/16. We have 
performed a numerical simulation for L = 200 under the no-flux boundary 
condition. The initial condition is u(x) = —l + 2x/L. Figure 2 displays the 
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Figure 2: Final stationary state for Eq. (5) at K\ = K 2 = 0.5, c\ = 
-0.1, c 2 = and L = 200. 
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final state in the time evolution by Eq. (5) at K\ = 0.5, K 2 = 0.5, c\ = —0.1 
and C2 = 0. There appear three domains of u ~ -x/3/2, and y/3/2. 
If the domain sizes are denoted as 1+,Iq and respectively for the three 
domains of u = v/3/2, and — y/3/2, the domain sizes satisfy approximately 
(u) = V3/2(l + -l_)/L = a = -0.1 and (u 2 ) = 3/4(Z + +Z_)/L = 3/{16K 2 )+ 
C2 = 3/8. The expected domain sizes are l + = 38.4, Iq = 100 and I- = 61.6 
for the parmeters c\ = — 0.1, C2 = and L = 200. The numerical result 
is approximately l + ~ 39.7, Iq ~ 97.5 and ~ 62.7. The global feedback 
succeeds in the domain size control for this quintic Ginzburg-Landau system. 

The same method is applicable for nonvariational systems. We use the 
quintic complex Ginzburg-Landau equation: 

dW n A d 2 W 

( ^- = - aW + {l + ic 2 )\W\ 2 W-\W\ A W + ^ r , (6) 

where W is a complex variable and C2 is a nonvariational parameter. There 
are two stable uniform states: W = 0, and W = Woexp(ia;ot) where Wq = 



1 + VI — 4a)/2 and loq = C2(l + \/l — 4a)/2. There exists a domain wall 
which connects the zero state and the oscillating state. Figure 3 displays 




Figure 3: Velocity v of the domain wall as a function of a for Eq. (6) at 
c 2 = 0.4. 



the numerically obtained velocity of the domain wall as a function of a for 
C2 = 0.4. The positive velocity implies that the oscillating state invades the 
zero state. The velocity of the domain wall takes at a = a c ~ 0.0678. The 
model equation with a global feedback is 

dW d 2 W 
— = -aW + (1 + ic 2 ) | W | 2 W- | W | 4 W - K(\ W \ 2 )W + (7) 
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where K denotes the feedback strength and (| W | 2 ) = (l/L) J L | W | 2 dx. 
Figure 4 displays the time evolution of | W | for a = 0.01, L = 400, c 2 = 0.4 
and K = 0.12. The initial condition is KeW(x) = 1 and Imlf(x) = for 
196 < x < 204, and W(x) = for x < 196 and x > 204. The domain size I 
of the oscillating state is calculated by the relation a + K(\ W | 2 ) ~ a + K \ 
W\ | 2 l/L ~ a c , where | W\ |~ 0.866 is the amplitude of the oscillating 
state coexisting with the zero state at a = a c for Eq. (6). The estimated 
value is I = L(a c — a)/(K \ W\ | 2 ) ~ 256.8. The numerically obtained 
size is approximately 257. We can control the domain size by changing the 
parameter a or K. 
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Figure 4: Time evolution of | W \ for Eq. (7) at a = 0.01, c 2 = 0.4, K = 0.12 
and L = 400. 



We have performed a simulation of the two-dimensional quintic complex 
Ginzburg-Landau equation with the global feedback. 

^ = -aW + (1 + ic 2 ) | W | 2 W- | W | 4 W - K(\ W \ 2 )W + V 2 W, (8) 

where K denotes the feedback strength and (| W | 2 ) = (l/L 2 ) Jq Jq \ W \ 2 
dxdy. The parameters are L = 200, a = —0.05, c 2 = 0.4 and K = 0.4. 
The numerical simulation was performed with the psudospectral method 
of mode-number 256 x 256 and timestep 0.005. The initial condition was 
ReW(x,y) = 0.0033(x - L/2)( r d -r) and lmW(x,y) = 0.0033(j/ - L/2)(r d - 
r) for r < rd where r = \J(x — L/2) 2 + (y — L/2) 2 and = 55, and 
W = for r > rd- That is, a topological defect is set at the center 
(x,y) = (L/2, L/2) as an initial condition. A spiral pattern evolves from 
the initial condition. The spiral pattern occupies only a finite domain by 
the global feedback. Figure 5 displays a 3D-plot of ReW(x, y) for the local- 
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Figure 5: 3D-plot of ReVF for a localized spiral pattern as a result of the 
time evolution by Eq. (8) at a = -0.05, c 2 = 0.4, K = 0.4 and L = 200. 



ized spiral pattern. The value of a + K(\ W | 2 ) is approximately -0.0648 for 
the stationary localized spiral, and the value is close to o c , but it is slightly 
larger than a c . It is probably due to the surface tension effect. Recently a lo- 
calized spiral pattern is studied with the quintic complex Ginzburg-Landau 
without the global feedback term [12]. However, the localization mechanism 



is different. In our model, the domain size of the spiral pattern can be 
controlled by the parameter a or K. 

To summarize, we have studied the control of the domain size for the 
Ginzburg-Landau type equations with a global negative feedback. The 
domain-size control by the global negative feedback is one of the simplest 
examples of the control for spatially extended dynamical system. The do- 
main size is determined by the relation that the velocity of the domain wall 
is zero. The global negative feedback can be derived from the coupled re- 
action diffusion equation. In the crystal growth problem, the temperature 
field plays a role of the inhibitory medium through the latent heat released 
at the solidification. The global negative feedback appears more naturally 



in electric circuits [13]. The control of the domain size is fairly robust and 



may be applicable for many systems. 
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